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Effect of Concentration on the Viscosity 
of Dilute Solutions* 

By Robert Simha 

A theory of the higher coefRcients of the viscosity-concentration curve of a suspension 
of nonspherical particles is presented. Starting with the simple model of a dumbbell, the 
flow around a single particle and its modification due to interaction with other particles are 
considered. It is shown that the coefficients a,- in the equation 

r} = 7jo{l + aiC+a2C^+azC^+, . . .) 

are related to each other, namely, 

a2=kiaij az=k2a^, . . . , 

where the /b's are independent of molecular size, in agreement with empirical equations. An 
explicit value for hi is obtained that applies also to models consisting of various arrays of 
spheres. The variation of hi with molecular shape is discussed, and various factors affecting 
the numerical values and the validity of the last-mentioned equation are pointed out. 



I. Introduction 

The dependence of the viscosity 77 of a dilute 
solution on the concentration c (expressed as 
volume fraction or as weight by volume) can be 
represented as a power series ot the type 



r7-=r7o(l + aiC+a2c2+a3c3 + ), 



(1) 



if the solute is a nonelectrolyte. 770 is the viscosity 
of the pure solvent. Considerable experimental 
and theoretical effort has been devoted to a 
determination of the coefficient ai, the so-called 
intrinsic viscosity, for suspensions of particles of 
various shapes since Einstein's [1] ^ original 
hydrodynamic treatment of spheres. In these 
hydrodynamic theories, ai is a measure of the 
disturbance of the flow of the solvent due to the 
presence of the solute particles at infinite dilution. 
The total effect under such conditions is the sum 
of the effects of the individual particles. For the 
range of higher concentrations, a number of 

1 Figures in brackets indicate the literature references at the end of this 
paper. 

♦Awarded an A. Cressy Morrison Prize in Natural Science in 1948 by the 
New York Academy of Sciences. 
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empirical equations have been proposed. They 
can all be expressed as a power series of type given 
in eq 1. As one can readily see, the theoretical 
calculation of the higher coefficients a2, <^3, etc., 
becomes increasingly complicated. It has been 
carried out exphcitly first for the coefficient as of 
a spherical suspension by the author [2]. ^2, a 
number independent of size for spheres, measiu-es 
the disturbance of the flow around a particle by 
a second one at a finite distance from the first. 
The method applied in the particular case is 
analogous to the procedure first developed by 
Smoluchowski [3] to obtain the modification of 
Stokes' law for the resistance of a sphere caused 
by the presence of a second one. Smoluchowski's 
results, oi which we sliaU make use, are obtained 
by a perturbation in terms of inverse powers of the 
distance between the two spheres. 

Empirically it has been suggested [4, 5] that the 
coefficients ai and a2 are related in polymer solu- 
tions by the equation 

a2=kai^, (2) 

where the factor k is independent of the dimensions 
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of the solute particle but may vary with the 
solute-solvent system. In several instances, 
although not always, this rule seems to apply 
satisfactorily. It is the purpose of this paper to 
examine the validity of the proportionality (eq 2), 
using first the simplest model for a nonspherical 
particle, namely, that of a dumbbell considered 
originally by Kuhn in his treatm.ent of intrinsic 
viscosity [6]. This model can be extended to 
that of a pearl string currently used by various 
authors in the description of hydrodynamic 
properties of chain molecules. 

In section II the velocity distribution around 
a single particle and the intrinsic viscosity ai are 
obtained because the results are needed for what 
follows. In sections III and IV, interaction 
effects between solute molecules are considered. 
In section V attention is devoted to the nature of 
the approximations used and to the effect of vari- 
ous factors on the validity of eq 2 and the magni- 
tude of the parameter k. 

II . Flow around an isolated particle. 
Intrinsic viscosity 

Figure 1 illustrates the notations and reference 
axes used. 2iis the length of the dumbbell and 
a the radius. It will be assumed that a/L<^l. 
Let the flow of the pure solvent be represented by 
a simple shearing motion in the XiX2-plane, 
namely, 



K=qX2, U2'=u,'=0, 



(3) 



where q is the velocity gradient and Ui^ the velocit}^ 
component parallel to the Xraxis. If the dumbbell 
is held in a given orientation relative to the refer- 
ence system Xi, the centers P and Q have the co- 
ordinates Xjc(P) and X}c(Q) = —Xk(P), respectively. 
Equation 3 then gives for the velocities at P and () 



u,{P) = qx2{P); u,{Q) = -u,{P), 



(3a) 



where the Xjc{P) are functions of the length L and 
the direction cosines of the particle axis. In the 
absence of rotatory Brownian motion, this flow 
will rotate the particle axis into preferred direc- 
tions. If one superimposes upon the motion (eq 3) 
a second one resulting from the spins of the par- 
ticle and considers the resulting force and torques 
acting on the spheres P and Q, one arriv^es at a 
result corresponding in three dimensions to 
Kuhn's [6] and identical with Jeffery's [7] formulae 




Figure 1. Dumbbell model and reference axes used and 
geometric relations involved. 

for the rate of change of the orientation specialized 
for a thin rod. Here however we shall be con- 
cerned only with the case of so-called complete 
Brownian motion, where any orientation tendency 
is absent. Following a previous procedure [8] we 
shall consider only the effects of the velocities 
(eq 3) and (eq 3a) on the particle. 

The change in the unperturbed flow (eq 3) due 
to a single dumbbell is obtained as a solution of the 
hydrodynamic equations of motion, which neglect- 
ing inertia terms, assume in this case the form 



'noV^Ui= 



^Un, 



^x„. 



FXi 
=0. 



(4) 



Here, as in the following, summations from one 
to three are performed over twice occuring indices 
m. ^ is the pressure. At infinite distance R, 
Ui has to assume the form shown in eq 3. If 
the velocity Ui is transformed to a reference sys- 
tem Xi{P) parallel to the first but with its origin 
at P, such that X2=X2(P) + X2(P), it is seen that 
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the velocities Ui consist of two parts if the mutual 
disturbance between the two spheres is neglected. 
The first satisfies the condition imposed by the 
X2(P) — term. The required solution is of the 
type given by Einstein [1]. The second part, 
satisfying the constant velocity qx2(P) is of the 
type given by Stokes [9] for translational flow. 
Now for the calculation of the dissipated energy 
of the system particle plus solvent, it is sufficient 
to consider the velocity and pressure distribution 
at large distances compared with the dimensions 
of the solute molecule, and therefore a fortiori 
large compared with the radii a of the spheres. 
As Einstein^s solution contains terms of the order 
a^/pp^ (multiplied by a velocity gradient) and 
higher, while in Stokes^ solution terms of the order 
a/pp (multiplied by a velocity) are found, the 
former can be neglected for our immediate pur- 
poses. Physically this amounts to concentrating 
the effect of each sphere on the solvent flow in its 
center P or Q. The calculation of the perturbing 
field of flow^ then bears a certain formal similarity 
to that of the field of a rigid dipole, and including 
hydrodynamic interaction, to the calculation of 
dipole-dipole interaction. It should be remem- 
bered, however, that at distances close to the 
particle the equations given below are invalid. 
For the disturbance due to sphere P, we find in 
this manner from Stokes' solution neglecting terms 
of the order a^/pp^ : 



„.-,-|„,OT°,,-|„,(P) «^.(f)f(^) 

4 Pp 4 Pp 

^=const. — 7T 7joUi{P)a ^\ j 

L Pp^ 



(5) 



where 






The additional velocity n? around sphere Q 
is given by an analogous equation. From eq 
3a and 5 it will be seen that the total velocity 
obtained by superposition of Ui^P) and UiiQ) 
will contain terms 



and 
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At points in the field of flow for which not only 
Ppj Po^dy but also R>L, it is permissible to 
expand these quantities in terms of the ratio 
L/R (see fig. 1) and to discard the third and higher 
powers of L. Omitting details, one finds 



112 
Pp Pq ri k 



Y (6) 



X,{P)X,{P) XmXi{Q) _ 

pp^ Pq' 

^-^IZ X,x,{P)-^,[X,x,{P) + X,x,{P)] 

PP^ pg3 - R^ Z^^uXAt") ^3 • 



The velocity to be added to eq 3 is then given by 

~~2 w y^'^^^k{P)^\i— 

[Jr.:r, (P) + X,x,(P)] + ^#'X;^,a;,(P) I ( (5a) 



u 



P=co„st.-3,.*(|J-«p2:j*m_„(P, 



Equation 5a is, of course, a solution of the equa- 
tions of motion (eq 4). It is useful to write it 
in the form of Lamb's [9] general solution in terms 
of spherical harmonics. It is easily found that 
eq 5a corresponds, using his notation, to: 



P-3 = P, x-2=~2 



--^^^PlXMP)-XMP)]- (50 



In order to obtain the dissipation of energy 
caused by the flow (eq 5a), one considers a large 

4:7r 

spherical volume V=^ Rq^ around the particle. 

Let prxi^ and Prx^ denote the components of the 
radial stresses exerted per unit area on the surface 
of this sphere. Then one derives from eq 3 and 
the definition of the stresses: 



V 

PrX,^=Voq~^'^PrX2^ 



--Voq^'-^Prx.'^O. (3b) 



Equation 5', together with an expression cited in 
[9] yield 

Prx = l^^o^^(-P)^ {g^^ J:X,x,{P)-XMP) } ' (7) 
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The rate of dissipation per unit volume is given by: 

= yj iPrxJ + VrxJ{uJ + Um)dS (8) 



dW 1 

dt 



where the integral is extended over the surface of 
the sphere Rq. The result is, omitting terms of 
the order 1/Rq and using eq.''3, 3a, 3b, 5a, and 7: 



dW 

dt 






if the particle is held in a fixed orientation. If the 
solution consists of n particles dissolved in the 
volume F, such that N=nlV represents their 
average number per unit volume, one has 



2^aW=c, 



(9) 



where c is the volume concentration. In the pres- 
ent case all effects are additive, and hence the 
dissipation of the whole solution is obtained by 
combining eq. 8a and 9 to 



dW 

dt 



-mMi 



(> 



27 X2{Py 
'10 a^ 



(8b) 



In order to obtain the viscosity r; of the solution 
we follow the method of reference [1].^ It will be 
assumed that the solute molecules are randomly 
distributed in the volume V. Let particle i have 
its center situated at a point ^k{i) in the X^-system, 
k=lj2,3',i=l,2, . . . 7^. Its contribution i/^C'i) 
to the flow disturbance at a point Xjc is obtained 
from eq 5a if the coordinates Xk are replaced by 
the relative coordinates Xk—^kii)- Hence we can 
write for the velocity distribution in the whole 
solution 



Ui=8uqX2+^Ui{i). 

i = l 



(5b) 



The viscosity of the pure solvent can be defined by 
the relation 

(,Y>*'1^^^*=''»2 -"oL^^+^J^^^; (10) 

In the same manner we regard the solution as a 
liquid with effective viscosity rj subject to an 
effective gradient g, defined as 



dW 

dt 



2 fbui , du2~\^ 



M^^t^m^'^lJ 



(10a) 



where the left-hand side is given by eq 8b. As 
d/dX;fc= — b/dffc(^), the sum in eq 10a may be 
replaced by 

JvV>^2^(>^Jxk=0 

if an integral is substituted for the sum. The 
integration is to be carried out over the whole 
volume V, The expression can be transformed 
into an integral over the surface S, namely, 



Substitution from eq 5a then leads to the result: 

J^^AilA (11) 

10 a' V ^ ' 



.(. 



Thus the viscosity tj of the solution follows from 
eq 10a as 

f- 27a:.(P)nr, 9 x,{Py Y 

Equation 12 gives the correct viscosity up to Unear 
terms in the concentration c. Averaging over-all 
(random) orientations, {x2{Py)={L^/S), we 
finally obtain 



77=770(1 + 2'^^)' 



(13) 



The intrinsic viscosity for the model under con- 
sideration is proportional to the ratio L^/a^ in 
concordance with the familiar result obtained for 
infinitely thin ellipsoids or rods not subject to 
any orientation effects by the flow gradient. 

The quadratic and higher terms in c obtained 
from eq 12 have to be discarded as long as hydro- 
dynamic interaction effects that give contributions 
of the same character and order of magnitude [2], 
are not taken into account. This has sometimes 
not been observed in the literature when Ein- 
stein's results have been extended to higher con- 
centrations. However, it can be seen that the 
higher terms in eq 12 will, after averaging the 
powers of X2(P), contain successively higher powers 



2 The procedure given in reference (7] can equally well be made use of when 
interaction effects are absent. 



412 



Journal of Research 



of the ratio i^/a^, thus confirming the empirical 
relation (eq 2) which may now be written as 



a2=kiai^; a^=k2ai^. . . ., 



(2a) 



where the fc's are, for the present model, pui'e 
nimibers and independent of the dimensions of 
the molecule. 

Equation 13 could, of course, have been derived 
in a much more straight-forward manner by con- 
sidering the frictional resistance to the flow 
offered by each of the two spheres. The rate of 
dissipation arising in this way is given by the 
product of frictional force and relative velocity. 
If the latter quantity is taken from eq 3a, we have 
by application of Stokes' law 



dt 



^m'+^^TVoaq" [x2{Py+X2(Qy] N, 



Elimination of A^ by means of eq 9 leads at once 
to eq 13. 

, III. Interaction between two particles 

The influence of a second particle on the viscous 
flow around a first has been considered first by 
Smoluchowski [3] for spheres in translational mo- 
tion, leading to a modification of Stokes' law. 
The effect so obtained is a function of the ratio 
between radii and mutual distance of the two 
spheres. More exact calculations are reviewed in 
Oseen's treatise [3]. Rotating spheres have also 
been treated [2] and calculations for nonspherical 
particles have been presented [3, 1 0]. The principle 
of Smoluchowski's method consists in seeking to 
fulfill by successive ''reflections" the hydrody- 
namic boundary conditions on the surfaces of all 
solid bodies involved. That is, additional flow 
distributions are superimposed on the previous 
flow, to compensate for the additional disturbance 
at the boundary between particle and fluid in- 
duced by the neighbors of the former. Thus the 
boundary conditions are fulfilled to successively 
higher degrees of approximation. The first ap- 
proximation that is sufficient for our present pur- 



pose has been shown [3] to fulfill the boundary 
conditions not at each point of the surface, which 
would make it an exact solution, but on the 
average over the whole surface. The velocity 
distribution is in other words such that, if aver- 
aged over the surface of the sphere, it gives no 
relative motion between the sphere and the liquid 
at their boundary. The solution previously em- 
ployed by the author [2] in the viscosity treatment 
is of the same degree of accuracy. 

Smoluchowski's original equations refer to 
spheres moving in the same direction with equal 
speeds. His results are modified slightly if this 
equality does not hold. Let the position of the 
centers of two spheres be given by ^^^^^ and ^^^^^ in 
a coordinate system fixed in space, their mutual 
distance by R^ and their respective velocities in 
the Ji-direction by Ui^^\ Ui^^K At a point ^^ 
(fig. 2) the velocity and pressure distribution pro- 




•*f2 



Figure 2. Geometric relations in interaction of two spheres, 
equations I4 and IS. 



duced by sphere 1 at large distances i?i has the 
following form : ^ 



j,(., = 3^^^(|™^™^ 



Bi' 



hWSm-U,'''U^,''Ki,'''] 



where 



(14) 



t^."'(?.)=f|«.+f|.ar 



■?;"') (?r 



■?/") 



i The author is indebted to E. O. Knox for checking eq 14. 
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■U ij^^^ {^]c^^^) represents this function at the point 



^.''\ 



An analogous expression results for Ui^^^ and p^^\ 
It is seen that Ui^^^ and p^^^ each consist of two 
terms. The first is, of course, Stokes' expression. 
The second one is a function of the -mutual distance 
and represents the interaction. Equation 14 
leads to a frictional resistance 



with 



[X^-x^{P)][X,-x,iP)]) 



pp 



\ 



(16) 



3 a 



(?i<^'-|i">) (?/''-^/") W'\ (15) 



Thus the resultant resistance is no longer in the 
direction of motion. 

Let us now consider a second dumbbell (M, ISl) 
as shown in figure 3. Relative to its center the 
points M and A^ have the coordinates X]c{M) and 
—Xjc(M), which are functions of the orientation of 
the axis MN in respect to the X^-system. They 
determine Ui{M) and Ui(N) — —Ui(M), the re- 
spective velocities of the undisturbed solvent flow 
relative to the particle at points M and A^. As in 
section II we shall be concerned with distances E, 
which are large in comparison with L. Consider- 
ing furthermore only molecules that are far apart, 
the modification of the Einstein-type of motion 
which gives rise to terms of the order a^ [2], will be 
omitted and only a^-terms retained. With these 
approximations one can make use of eq 14 and 
write the velocity around P to be superimposed 
on u/' in eq 5 



Ui/(M, P) is obtained from the last equation by 
replacing pp by ppm and the coordinates Xj by 
^j-\-Xj{M). The last term in eq 16 is obtained in a 
corresponding manner. An expression analogous 
to eq 16 results for Aui^. We define the inter- 
action coefficients pji and g^i by means of the 
relations: 

J),, = Uj,^iM, P) - ?7,/(iV, P); q,, = L\<^{M, Q) - 

Un''{N,Q) (17) 

and obtain from eq 16 for the total velocity dis-' 
tribution around the molecule P, Q 

ui=ui (Stokes) +u,{M)pmiUiM''{X,P) + 



U,{M)qr,aUi.r.%X,Q) 



(16a) 



where the first term is given by eq 5a. The re- 
mainder contains functions of the mutual dis- 
tance and orientation of the two interacting mole- 
cules. As stated previously, one sums over m 
from one to three. 



The quantities Uim^ and Uim^ in eq 16a are of the same type as Stokes' velocity distribution and can 
be transformed in the same manner for large distances E. One finds in this manner in analogy to eq 5a 



'Z/i='Ui(Stokes)-t-^ — ^3 — {jpm\ — q_m\) 



E^ 



^=i> (Stokes) +2^0 — ^3— 



X,a;,(P)5,,-[X^xKP) + X.:^.(P)] + ^^^Z;-^.^^(^) 



{P.nl-qr.l) {^-§'i:X,X,{P)-X^(P)y 



(16b) 



In analogy to eq 7 this makes a contribution to the J^i — component of the radial stress 
i?rxi = i>rxi (Stokes)- 2 r;oa — gr" {Vmi—q,mi){- 



E' 



^ S XjMP)-[XnMP)+X^Xm{P)] 



(18) 



Application of eq 8 now yields for the total rate of dissipation per unit volume of the system consisting 
of solvent plus one particle (P, Q) acted upon by a second one (M, A^): 



dW 

dt ~ 



= M 



1 +f y X.iPy-l'^ x,{M) [3X2(P) {pn -qn)-2x,{P) (P2i-g2i)]|- (19) 
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Figure 3. Geometric relations in interaction between two 
molecules J equation 16. 

The interaction coefficients p^i and g^i can be 
evaluated for pairs sufficiently far apart in terms 
of the mutual distance p (see fig. 3). A series 
expansion in terms of 1/p gives, using eq 17 and 
discarding higher powers in i/p: 

x,{P)xj{M) + x,{P)x,{M) + bn [^ Z:te(P). 

^i,x,{M) - j:x,{P)x,{M)j j . (20) 

Equations 16b and 20 constitute the complete 
solution of the interaction problem for two 
particles. 

IV. Viscosity of Solution 

To obtain the quadratic concentration term in 
the viscosity one computes first from eq 19 the 
dissipation of the mixture. If we visualize again 
the centers of the solute molecules randomly dis- 
tributed over the volume Vj then for the ith 
particle (PiQi) situated at a point ^kii) and acted 



upon by a fixed particle {M, N), p ji{i) — ^n{i) 
becomes a function of the relative coordinates 
ik—ik(S)' It is to be summed over all (PiQi) 
species interacting with one (M, N) and then mul- 
tiplied by the total number of particles. In this 
manner eq 19 yields 



dW 

dt 



JdW\ 

\ dt /stokes 



: rioTC 



X2{M) 



20 '"" a 

2Xl(P)Sb2l(^)-g2l(^)]}, 



^X2{Pm[Pn{i)-^n{i)]- 

(19a) 



where the first term is given by relation (eq 8b). 
The summations may be replaced by integrations 
over the volume V, The use of the expression 
eq 20 in the whole region amounts to neglecting 
the influence of particles close to a given one. 
The error so committed increases, of course, at 
higher concentrations where the average distances 
between solute molecules become small. The two 
integrals can easily be transformed into surface 
integrals by means of Gauss' theorem. Omitting 
all terms that vanish on integration because of 
reasons of symmetry, one finds for instance, that 
the remainder in (^21— <Z2i) ^^^ be written as 

(:P2i-g2i)--Y[:ri(P)x2(M) + a:2(P)xi(M)] 

with a similar expression for Pn—qn. The final 
result for the dissipation due to particles oriented 
in a given direction Xjc(P) interacting with particles 
in another direction Xk(M) is 



dW 

dt 



27 



l+?I^^:P-%- 



10 a" 
50^4 [^x,{P)x,{P)x,{M)x,{M) + 

QX2{P)X2(M) {X2{P)X2{M) + 

Xs{P)xs{M)) + x,(P)%{My] c' 



(19b) 



The next step consists in calculating the effec- 
tive velocity gradient g^ defined in eq 10a. The 
velocity distribution around molecule i in the 
presence of a perturbing molecule at an arbitrary 
point ^jc is derived from eq 16b and 20 by re- 
placing Xjc by Xjc—^k{i) and ^jc by ^k—^kii)- 



Viscosity of Dilute Solutions 



415 



The corresponding result for the whole mixture 
follows again from eq 5b. g^ is then derived as in 
section II. The expression 



Ui r>-+t^2 



?1 



p r M/2 p 

Uq JXo 

to be integrated over the sphere Rq now involves 

fc 6 t 4 

terms of the order 4^ and wg- The integration 

JXq Kq 

requires the values of the following integrals over 
the surface of the spherical enclosure: 



fh%'dS=^^ R,^; f^,%HS^ 



%R,'; fii%%'dS=^^Rl 



The integrand written above depends upon the 
mutual orientation in each molecule pair. On 
averaging independently the X]c{P) and Xjc{M) 
terms, there remains after a lengthy but straight 
forward calculation, omitting all terms which 
vanish for random orientation: 

Instead of eq 12, one finds now, using eq 19b, tliat 

r 27 Xjjpy _27 7 x2{P)%{Mycn 

L "^10 a' ^ 10 ■ 5 a* J 



V = V(> 



[ 



, . 9 x,{Py 

10 a^ 10 



9 27 X2{Pyx2{My 



'1 



(12a) 



where the denominator represents qj^. Expand- 
ing up to quadratic terms in the concentration 
and averaging over the direction cosines, we arrive 
at the following result for the viscosity of the 
solution: 

On the other hand, eq 12 would have given for 
the coeflftcient of the quadratic term the smaller 
value 81/100. 



V- Discussion 

Equation 13a yields for the coefficient ki defined 
in 2a the value 

^1 = 0.77. (21) 

Had interaction effects not been taken into 
consideration, eq 12 would have given 



ki = 0.36. 



(21a) 



These results refer so far to the model of a dumb- 
bell. Let the molecule be represented next by a 
linear array of length 2i, consisting of a number 
of spheres with radii a, v per unit length, arranged 
evenly. The simple considerations at the con- 
clusion of section II, which lead to eq 13, indicate 
that now the factor L ^ is replaced by an average 
value {U)<IU. This is seen as follows: The 
contribution to the dissipated energy by a sphere 
located at a distance y from the center of the 
molecule is proportional to y^. This is to be 
summed over all groups in the molecule and 
multiplied by N^ the number of molecules in unit 
volume. A^ is inversely proportional to 2vL, 
Thus the intrinsic viscosity ai becomes propor- 
tional to 



^J>y=(i:^). 



In the case of a randomly coiled chain one aver- 
ages over all internal configurations of the chain in 
addition to averaging over all orientations of the 
whole molecule. As long as coiling is not extensive 
the determining quantity has been shown re- 
peatedly to be the mean square separation of 
chain ends. Reverting to the linear array, it can 
be seen from eq 19 and 20 that the interaction of 
parts of two molecules at distances yi and ?/2 from 
their respective centers will contribute terms pro- 
portional to yi^y2^ to the dissipation rate. Con- 
sequently, the coefficient ^2 will be proportional to 
the square of {U) as one could have surmised. 
A corresponding result obtains for a solution con- 
sisting of moderately coiled chains. Thus the 
numerical values (eq 21 and 21a) should be valid 
also for a long extended molecule and a coil. 

The value (eq 21) is larger than those derived 
experimentally from suitable extrapolations of 
viscosity-concentration curves, the latter being 
usually of the order of 0.3-0.7 and more nearly 
comparable with the value (eq 21a). On the 
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other hand, the theory in its present form can 
reproduce only approximately the contribution of 
the hydrodynamic interaction to the ki coefficient 
or rather to its excess over 0.36. It has been 
stressed previously that the evaluation of the 
interaction coefficients pji, qji (eq 17) as a series 
(eq 20) and therefore eq 19b are valid only for 
large intermolecular distances. One could im- 
prove the situation by dividing the volume V 
around the center of a molecule into three parts. 
The first comprises all particles with spherical 
groups so close to those of the center particle, that 
an expansion of the mutual disturbance in terms 
of the ratio alRi2 (fig. 2) becomes altogether 
meaningless. In other words, solute molecules are 
practically touching each other whatever their 
relative orientation. From a hydrodynamic point 
of view, these complexes would be treated as 
single molecules. Such conditions, frequent at 
higher concentrations for large molecules will be 
disregarded here (see also below). In the second 
region intermolecular distances are large in com- 
parison with the size of the individual group in the 
molecule, but are of the order of magnitude of the 
dimensions of the molecule. This would permit 
the application of Smoluchowski^s equation (eq 
14), but not of the expansion (eq 20). All 
further summations are to be performed in terms 
of individual group to group distances without 
using expansions in terms of molecular center to 
center distances. In the third region the current 
approximations are valid. Of course, the sub- 
division is somewhat arbitrary. The conditions 
in region two affect both the numerator and 
denominator in eq 12a. They should be more 
important in the calculation of the dissipated 
energy (see eq 19a) than in that of the effective 
velocity gradient qs (compare eq 10a), because of 
the direct summations over the {Vn~~^n) occur- 
ring in the former. Also in this region a disturb- 
ance of the random mutual orientation of particles 
merely because of geometric reasons is possible. 

The treatment in section II took no account of 
the hydrodynamic interaction between groups 
in the same molecule. For the dumbbell a simple 
calculation based on the result (eq 15) leads in a 
first approximation to a correction factor in the 
expression (eq 13) for Oi equal to 



1 + 



20 L 



for ajL <C 1 . Also as a first approximation , ^2 remains 
unaltered. Thus only a decrease of ki by a few 
percent results from this effect. In an array of 
spheres the interaction of neighboring spheres 
v\dth velocities of the same sign relative to the 
surrounding liquid will actually give a correction 
factor that is smaller than unity. The inter- 
action between different groups within a strongly 
coiled molecule changes the flow pattern pro- 
nouncedly, as has been shown recently by several 
authors, and this influences also the total mutual 
disturbance between two different molecules. 

It is perhaps noteworthy that the calculation in 
section II already gives as large a value of ki as 
0.36. Effects other than purely hydrodynamic 
ones may contribute to the slope of the viscosity- 
concentration curve even at relatively moderate 
concentrations. Experimentally the reported ini- 
tial slopes have sometimes been obtained by extra- 
polation at concentrations at which the conditions 
inherent in our derivation are no longer fulfilled; 
the solute molecules of large molecular weights 
are no longer sufficiently far apart on the average. 
There are indications in the literature that at 
concentrations roughly below 0.1 g/100 ml solu- 
tions of chain polymers exhibit steeper slopes than 
the extrapolated straight lines. Qualitatively 
similar effects have been observed in tobacco 
mosaic virus solutions [11]. In the latter case one 
is dealing with long rod-like particles and the 
possibility of aggregation is strongly increased at 
higher concentrations. Orientation effects are then 
relatively more important and this reduces the 
viscosity at higher concentrations. Systematic 
measurements of the viscosity of very dilute 
polymer solutions as function of velocity gradient 
and solvent character are desirable to decide the 
nature of the apparent downward curvature, 
found on approaching the region of high dilution. 

It is obvious that the present considerations do 
not apply directly if the nature or the identity of 
the individual molecule changes with concentra- 
tion. For instance if formation of aggregates 
occurs which can be expressed, say, by an equi- 
librium constant K between double and single 
molecules it is seen that this m.akes an 
additional contribution to the coefficient a2, namely, 

where ai^^*— ai is the difference between the 
intrinsic viscosity of double and single molecules. 
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Thus the higher coefRcients in the expansion 
(eq 1) cannot be found solely by hydrodynamic 
considerations and therefore eq 2a is then not 
valid. This is particularly the case when the 
effect of velocity gradient in concentrated solutions 
is to be examined. However, even in dilute 
solution the interaction is more affected by 
orientation than the intrinsic viscosity. 

Finally the present results for asymmetric 
particles may be compared with those previously 
obtained for compact spheres. The author's 
method of calculation [2] leads to a value (i2=14.1 
[10]. Einstein's results give a value of 4 for a2. 
Thus, for spherical particles, the corresponding 
values are /:i = 2.26 and 0.64, respectively. It 
would be of interest to calculate the gradual 
increase in the Z^i-values with increasing sym- 
metry by considering ellipsoidal particles of finite 
thickness. Such results would be of value in the 
case of protein solutions and also of chain poly- 
mers. One should expect a parallelism between 
ki and other parameters characterizing shape 
such as the exponent in the modified Staudinger 
equation. Also branching and cross-linking 
should result in an increase of ki. Such changes 
have been reported. It appears from the fore- 
going considerations that, leaving aside other fac- 
tors, the parameter ki will remain independent of 
molecular weight as long as the degree of coiling 
does not change markedly. 

VI. Conclusions 



The numerical value derived for ki (eq 2a and 21) 
is of the correct order of magnitude but somewhat 
larger than those encountered for solutions of chain 
molecules in good solvents, to which the model 
used is applicable under certain conditions. It 
seems to be more nearly in agreement with some 
numerical results for globular proteins. Possible 
experimental and theoretical reasons for this 
divergence have been outlined. 

The parameter ki depends in a characteristic 
manner on the shape, as determined for coiling 
molecules by their structure, the solvent environ- 
ment, and also on the molecular weight. This 
conclusion can be reached on the basis of the 
results obtained for two extreme cases treated here 
and in previous work, respectively. In this con- 
nection it is noteworthy that globular proteins 
appear to give larger values for ki than chain 
polymers. The region between the two extremes 
of an '^open'' coil, or dumbbell, and a sphere re- 
mains to be treated. On the basis of such a 
treatment it may become possible to derive 
information about the properties of the solute 
similar to that presently obtained from intrinsic 
viscosities. 

It is of interest to examine the influence of the 
mutual disturbance on the viscosity at finite 
velocity gradients, when partial orientation in the 
field of flow occurs. Effects present at infinite 
dilution should be more pronounced. 

VII. References 



In treating the concentration dependence of 
r]sp, two effects may be considered at low concen- 
trations. Aggregation, or entanglements, and 
hydrodynamic interaction of distinct molecular 
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viscosity of the solution, which depends on the 
frequency of the aggregates and can be described 
in terms of the intrinsic viscosity of an aggregate. 
The second one leads to results that are in agree- 
ment with empirical relations proposed for dilute 
solutions. Without explicit calculation of the 
cubic and higher terms in c it cannot be seen 
whether -qsp is represented by the binomial for- 
mulas given by various authors. The present work 
indicates that they probably are not strictly valid. 
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